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ABSTRACT 

We describe the dynamical evolution of hot spots velocity, pressure and mass density 
in radio loud active galactic nuclei (AGNs), taking proper account of (1) the con- 
servations of the mass, momentum, and kinetic energy flux of the unshocked jet, (2) 
the deceleration process of the jet by shocks, and (3) the cocoon expansion without 
assuming the constant aspect ratio of the cocoon. By the detailed comparison with 
two dimensional relativistic hydrodynamic simulations, we show that our model well 
reproduces the whole evolution of relativistic jets. Our model can explain also the 
observational trends of the velocity, the pressure, the size, and mass density of hot 
spots in compact symmetric objects (CSOs) and FR II radio galaxies. 
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1 INTRODUCTION 

Which evolutionary tracks are radio loud AGNs (radio 
galaxies) passing through? This is one of the primal issues 
in the study of AGNs (Ryle & Longair 1967; Carvalho 1985; 
Fanti et al. 1995; De Young 1997). Stimulated by the ob- 
servational progress (e.g., Turland 1975; Readhead, Cohen 
& Blandford 1978; Bridle & Perley 1984), a number of hy- 
drodynamic simulations of jet propagations have been per- 
formed to examine their physical state of the jet (e.g., Nor- 
man et al. 1982; Wilson & Scheuer 1983; Smith et al. 1985; 
Clarke, Norman & Burns 1986; Lind et al. 1989; Clarke, 
Harris & Carilli 1997; Marti et al. 1997). These numerical 
studies have confirmed that the jet is composed of "light" 
(i.e., lower mass density) materials compared with an am- 
bient medium to reproduce the observed morphology of the 
expanding cocoon (e.g., Norman et al. 1982). However it is 
hard to examine the whole duration of powerful radio loud 
AGNs with sufficiently large dynamical range because of the 
limitation of computational powers. 

A new population of radio sources so-called "compact 
symmetric objects (CSOs)" has been recently noticed. The 
CSOs was first identified by Philips & Mutel (1980, 1982) 
and more complete sample were presented by Wilkinson et 
al. (1994) and Readhead et al. (1996a, b). Concerning the 
origin of CSOs, two scenarios were initially proposed. One 
is so-called "frustrated jet scenario" in which the ambient 
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medium is so dense that jet cannot break its way through, so 
sources are old and confined (van Breugel, Miley & Heckman 
1984). The other is "youth radio source scenario" in which 
CSOs are the young progenitor of FR II radio galaxies (e.g., 
Shklovsky 1965; Philips & Mutel 1982; Carvalho 1985; Fanti 
et al. 1995; Begelman 1996; Readhead et al. 1996a; O'Dea 
& Baum 1997). Recent observations reveal that their speeds 
are better understood within youth radio source scenario 
because of their age with 10 3_5 yr, which is much shorter 
than the age of FR II sources with 10 6_7 yr (e.g., Owaiank, 
Conway & Polatidis 1999; Murgia et al. 1999; Taylor et al. 
2000). This indicates the possibility of CSOs as the progen- 
itor of FR II sources although their evolutionary tracks are 
poorly understood. 

The hot spot, which is identified as the reverse shocked 
region of the decelerating jet, is one of the most important 
ingredients in the whole jet system. The evolution of the hot 
spot is tightly linked to that of cocoon because the cocoon is 
consist of the shocked plasma escaped from the hot spot (see 
Fig. 1). Observationally, the correlations between the hot 
spot properties (the velocity, the pressure, the size and the 
mass density) and projected linear size have been reported 
for CSOs and FR II sources (Readhead et al. 1996a; Jeyaku- 
mar & Saikia 2000; Perucho & Marti 2002 ). These obser- 
vational trends would also reflect the evolutionary tracks of 
radio loud AGNs. Thus, in order to clarify the physical re- 
lation between CSOs and FR II sources, it is inevitable to 
model the dynamical evolution of hot spots in radio loud 
AGNs. However, little attention has been paid to this point 
in spite of lots of theoretical evolutionary models have been 
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proposed based on cocoon dynamics (e.g., Falle 1991; Begel- 
man 1996; Kaiser & Alexander 1997). Thus, the goal of this 
paper is to construct an appropriate dynamical model of hot 
spots in the radio loud AGNs. 

The plan of the paper is organized as follows. In §2 
we outline and model a dynamical evolution of hot spots 
connected with the cocoon dynamics. In §3, we compare with 
previous theoretical and observational works. Conclusions 
are given in §4. 



2 DYNAMICAL EVOLUTION OF HOT SPOTS 
CONNECTED WITH THE COCOON 
EXPANSION 

2.1 Outline 

In this paper, we model a dynamical evolution of hot spots 
with the aid of cocoon dynamics (Begelman & Cioffi 1989: 
hereafter BC89; Kino & Kawakatu 2005: hereafter KK05). 
Specifically, the evolution of the hot spot velocity («hs), the 
hot spot pressure (Pus) and the hot spot density (phs) are 
discussed. These quantities are described in terms of the 
length from the center of the galaxy to the hot spot (Zh)- 

Concerning vns, radio observations of powerful FR II 
radio galaxies show us that hot spots are always reside at the 
tip of the radio lobe (e.g., Myers & Spangler 1985; Readhead 
et al. 1996b). Thus, it is natural to impose the relation of 



vhs = «h, 



(1) 



where Vh is the advance speed of the cocoon head. The veloc- 
ity Vh is significantly affected by the two dimensional (2D) 
effect. However it can be reasonably handled by the phe- 
nomenological description as follows. Consider a pair of jets 
propagating in an ambient medium (see Fig. 1). At the hot 
spot, the flow of the shocked matter is spread out by the 
oblique shocks that then deflects (Lind et al. 1989), the vor- 
tex occurs via shocks (e.g., Smith et al. 1985) and/or the 
effect of jittering of the jet (e.g., Williams & Gull 1985; Cox 
ct al. 1991) which behaves like the "dentist drill" (Scheuer 
1982). Thus, the effective "working surface" for the advanc- 
ing jet is larger than the cross section area of the hot spot 
Aj, which was pointed out by BC89. BC89 introduced the 
effective cross section area of the cocoon head Ah as that of 
the effective "working surface". Thus, we can determine 
the reasonable value of v b by the expanding cocoon process. 

As for Phs and pHS, we deal with them through one di- 
mensional (ID) shock junctions. Since the hot spot is identi- 
fied with the reverse shocked region of the jet, Phs and phs 
can be obtained as a function of Vh by combining with eq. 
(1)- 



2.2 Un-shocked Jet 

Here we set up the mass, momentum and energy flux of the 
jet with three assumptions. Our main assumptions are as 
follows; (i) We assume that the speed of the jet is relativis- 
tic on the large scale (~ 100 kpc) and the jet is consist of 

1 Before BC89, the head advance velocity is estimated by purely 
the ID momentum balance (e.g., Begelman, Blandford, and Rees 
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Figure 1. A schematic picture of the co-evolution of hot spots 
and cocoons. Most of the kinetic energy jet is injected via the 
termination shock of the jet which is identified as the hot spot. 
The sideways expansion speed of cocoon is v c . The area of the 
radio lobe at the position of hot spots is larger than that of 
hot spots. The head part of the cocoon advances with speed v^. 



the cold medium. Although the jet speed on large scales is 
still open issue, several jets are suggested to be relativistic 
ones (e.g., Tavecchio et al. 2000; Celotti, Ghisellini & Chi- 
aberge 2001; Uchiyama et al. 2005), (ii) The mass, energy 
and momentum of jets are conserved in time. Namely, we do 
not include the entrainment effect of the ambient medium. 
This is justified by the numerical simulations for highly rel- 
ativistic jet flows (e.g., Scheck et al. 2002; Mizuta, Yamada 
& Takabe 2004) and (iii) We ignore the dynamical effect of 
magnetic fields because of the kinetic flux dominance in FR 
II radio galaxies (e.g., Hardcastle & Worrall 2000; Leahy & 
Giani 2001; Isobe et al. 2002). Then, the mass (Jid), energy 
(Lid) and momentum (Qid) flux are given (Blandford & 
Rees 1974); 



Jid = TjAjPjC, 
Lid = r^AjpjC 3 , 
Qid = rjUjpjc 2 , 



(2) 
(3) 
(4) 



where Tj and pj are the Lorentz factor and the mass density 
of the jet, respectively. Note that the kinetic energy flux 
Lid denoted here satisfies the relation of Lj = (A h /Aj)LiD, 
where Lj is the total kinetic power shown in KK05. 

From these conditions, the following quantities are con- 
served for any l h ; 



Tj = const, 



PjAj c = const. 



(5) 
(6) 



1984). 



Interestingly, the Lorentz factor Fj does not depend on l^. 
In other words, the speed of jet is relativistic even on the 
large scale. 



2.3 Shock junctions between the jet and ambient 
medium 

We briefly review the ID shock jump conditions which gov- 
erns the deceleration of the relativistic jet by the surround- 
ing ambient medium (Kino & Takahara 2004 for details). 
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We can determine P a and p a from X-ray observations where 
P a and p a are the pressure and the mass density of the am- 
bient medium, respectively. The assumption of a cold jet is 
written as Pj = 0. We regard Fj as a parameter. 

Since vhs is estimated to be in the range 0.01c to 0.1c 
both for FR II sources (Liu, Pooley & Riley 1992; Scheuer 
1995) and CSOs (Conway 2002 and references therein), the 
forward shocked (FS) region quantities are determined by 
the shock jump conditions in non-relativistic limit (Lan- 
dau & Lifshitz 1959). By using the pressure balance along 
the contact discontinuity between the hot spot and ambient 
medium, we can obtain the expression of Phs as functions 
of two observable quantities vhs and p a such as 



Phs 



[5 - (l/M 2 )} , 
tP>'H<- 



15 [1-{1/M 2 )] 2 



(7) 



where M = vfs /^(5P a /3p a ) and vfs are the Mach num- 
ber and the velocity of the upstream of FS, respectively. 
We adopted the adiabatic index of the downstream of FS 
as 5/3. In the reverse shocked (RS) region, we employ the 
relativistic shock jump conditions in the strong shock limit 
(Blandford & McKee 1976). Then, the equation of state and 
the mass continuity in the RS region can be written as 



Phs = 



Pi = 



3Phs 



(Fj 



- l)c 2 ' 
3P HS 



(4^+3)^-1)^' 



(8) 
(9) 



where we set the adiabatic index in the RS region as 4/3. 
Thus, pHS and pj also can be given by p a and vhs- 



2.4 Dynamical evolution of the cocoon 

To determine the velocity of the cocoon head Vh with con- 
sidering 2D sideways expansion, we prepare the solutions 
of cocoon dynamics based on KK05. In KK05, by solv- 
ing the equation of motion along the jet axis, perpendic- 
ular to the axis (i.e., sideways expansion), and energy in- 
jection into the cocoon, we obtained the v c , Vh, P c , and 
in terms of ih, where v c , and P c are the velocity of co- 
coon sideways expansion and the pressure of cocoon, respec- 
tively. The declining mass density of the ambient medium 
is assumed to be p a (d) = p a o(d/do)~ Q , where d, do and 
pao are the radial distance from the center of the galaxy, 
the reference position and the mass density of the ambient 
medium at do, respectively. In order to convert i-dependence 
of the results of KK05 to independence, we use the equa- 
tion lh — J Q Vh(t')dt' and l c = J Q v c (t')dt' where l c is the 
radius of the cocoon body. The obtained cocoon quantities 
in KK05 are as follows; 



Vc = Vc0 {hTo) 

Pc=Pco(» 
VthO/ 



-2 + 0.5a) + 3 



X(l-q/2)-2 
l h \ X(-2 + 0.5a) + 3 



2-X(2-0.5a) 
X(-2 + 0.5a) + 3 



Vh = VhO 



(-) 

Wo J 



X(a-2)(-2 + 0.5a) + 3o 
X(-2 + 0.5q) + 3 



(10) 

(11) 

(12) 
(13) 



where X is the power law index of the effective cross section 
area of the cocoon body A c (t) oc t x (see Fig.l). Throughout 
this paper, we use the normalization of l b0 = «ho^7, where 
tj = 10 7 yr and Who is the velocity of the cocoon head at Iho 
and the coefficients of each physical quantities are denoted 
with the subscript 0. In KK05, we selected A c as an unknown 
parameter because we could not obtain the solution for a = 
2. However, by compare with previous works, it is worth to 
show the power law index f3 of the effective cross section 
area of the cocoon head A^ = Ahoilh/lho) 13 as a function of 
q and X. From eq. (13), we obtain the following relation as 



X(a - 2)(-2 + 0.5a) + 3a - 4 



(14) 



X(-2 + 0.5a)+3 

From that, it is clear that the head growth can be expressed 
accurately in KK05 even for a = 2. 

Self-similar models are also useful tools to explore the 
evolution of the cocoon expanding (e.g., Falle 1991; Begel- 
man 1996; Kaiser & Alexander 1997). However, the problem 
has been pointed out on the assumption of the constant as- 
pect ratio of cocoon employed in self-similar models. By the 
comparison of the young and grown-up sources, they claim 
that the condition of the constant aspect ratio is not ful- 
filled incidentally (e.g, Readhead et al. 1996a; De Young 
1997; Komissarov & Falle 1998; O'Dea 1998; Scheck et al. 
2002; Carvalho & O'Dea 2002a, b; Tinti & De Zotti 2006). 

2.5 Dynamical evolution of hot spots 

Combining with (i) the mass, momentum, and kinetic energy 
of the jet (§2.2), (ii) the deceleration process of the jet by 
shock (§2.3), and (hi) the cocoon expansion (§2.4), we can 
finally obtain the dynamical evolution of hot spots. From 
cqs. (1), (5), (7), (8), (9) and (12), the quantities vhs, Phs, 
Phs, and pj as follows; 



«hs = vhso 



(-) 

(^) 



2-X(2-0.5e») 
X(-2 + 0.5a) + 3 



Phs = Phso 



X(2-0.5q)(q-2) + 4-3q 
X(-2 + 0.5a) + 3 



X(2-0.5q)(ci-2)+4-3o 
X(-2 + 0.5a) + 3 



Pi = Pio 



(-) 

Wo J 



X(2-0.5a)(a-2) + 4- 
X(-2 + 0.5a) + 3 



(15) 



(16) 



(17) 



(18) 



Note that Phso, Phso and pjo can be expressed by only ob- 
servable quantities p a o and whso if 

we assume Fj (see eqs. 
(7), (8) and (9)). Thus, it is possible to know not only 1^ de- 
pendence but also the absolute quantities of hot spots and 
jets. The aspect ratio of the cocoon 1Z = l c /lh is the in- 
triguing quantity for studying the dynamical evolution of 
hot spots. The Zh-dependence of the aspect ratio of cocoon 
is then given by 



TZ = TZo 



Vtho/ 



X(2.5-0.5a)- 



2 + 0.5a) + 3 



(19) 



where TZo = (v c o/vhso)[(X(-2 + 0.5a) + 3)/(0.5X)]. As a 
consistency check of our assumption of constant A h /Aj, we 
can easily find that hot spot radius ths (oc A^ 2 ) shows the 
same Zh-dependence as A^ 2 from eqs. (6) and (18). From 
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above results, we obtain the slope of all physical quanti- 
ties as functions of two key physical quantities, namely a 
(the slope index of the ambient matter density) and X (the 
growth rate of cross section of cocoon body) . In the case of 
constant 1Z, our results agree with self-similar models of co- 
coon expansions (e.g., Falle 1991; Begelman 1996; Kaiser & 
Alexander 1997). However, we stress that these self-similar 
models assume that Phs/Pc and 7Z are both constant in lh, 
whilst we do not impose these assumptions and also predict 
the dynamical evolution of phs and pj. 

The relation between Phs and P c is also the interesting 
topic. From eq. iJJJ , the hot spot pressure is written by Phs = 
4/3a(^h)«Hs/3 for M 3> 1, while the over-pressured cocoon 
requires P c = p a (Z c )i>c- Thus, the ratio of Phs to P c is 



Table 1. Comparison with 2D hydrodynamic simulations and 
self-similar models 



Pc 

Phs 



0.5X 



X(-2 + 0.5a) 



(20) 



This implies that P c /Phs is controlled by 1Z and a. Since 
1Z < 1 and TZo < 1 are satisfied by definitions, P c /Phs 
should less than unity for < a < 2. In the case of 1Z = 
TZo = const or a = 2, it reduces to the interesting relation 
of 

Pc/Phs = 

This shows that P c /Phs is determined only by TZq. We stress 
that our model predict that P c is smaller than Phs as long 
as TZo < 1. Additionally, rewriting of the explicit form of 
the P c in terms of the quantities of the jet may be also 
stimulating, which is given by P c = TZoT?pjC 2 for IZo = 
const. From this, one can find that the larger Fj leads to the 
larger P c which is actually seen in relativistic hydrodynamic 
simulations (see Fig. 5 in Marti et al. 1997). To comprehend 
the energy injection process into the cocoon via the hot spot 
with the duration of ti n j, we rewrite the eq. 11201 as 



PcO^'cO'S'cO ~ PhSO^'cscO'S'hSO = Ljtin 



(21) 



where Shso = 47rr- HS0 , and S c o = 2irl c0 lho, ?Wo = 
c/(2(0.5X) 2 ) ~ (0.5 - 0.7)c. This describes the continu- 
ous energy injection of AGN jets (i.e., tinj = tage). On the 
contrary, Blandford and Rees (1974) used the relation of 
Pco(f c oS'co) 7c ~ -Phso^cscoShso) 7 -^ where 7 is the adia- 
batic index in each region. We claim that this relation is 
appropriate for the instantaneous (i.e., t[ n j <C ta ge ) injection 
seen in supernovae (SNe) or gamma-ray bursts (GRBs). 



3 COMPARISON WITH PREVIOUS WORKS 

3.1 Comparison with numerical simulations 

Scheck et al. (2002; hereafter S02) examined the long term 
evolution of the powerful jet propagating into a uniform am- 
bient medium (a = 0) with "2D" relativistic hydrodynamic 
simulations. S02 showed that the evolution of the jet pro- 
ceeds in two different phases appear (they are shown in Ta- 
ble 4 and Fig. 6 in S02). "ID" phase: In the initial phase 
(t < 1.2 x 10 5 yr), the jet shows ballistic propagation with 
Ah = const and whs = const. "2D" phase: This phase starts 
when the first large vortices are produced near the tip of 
the jet. Then, the cross section area of the cocoon head be- 
gins to increase. The hot spot then starts decelerating, but 
the jet speed remains the same relativistic one during whole 






"HS 


A h 


Pc 




Pi 


n 


"ID" Phase" 
S02 


.—0.11 


const 


7—0.95 


const 


const 


7-0.45 

'h 


BC89 


const 


const 




const 




,-0.5 

'h 


This work 


const 


const 




const 


const 


7-0.5 
'h 


"2D" Phase 6 
S02 


.-0.55 
l h 


;0.90 
l h 


7-1.30 
l h 


l h 


,-1.0 
'h 


7-0.09 
l h 


B96 


.-2/3 
'h 


.4/3 
'h 


,-4/3 
l h 


,-4/3 
'h 




const 


This work 


,-0.56 
'h 


I 11 

'h 


7-1.30 
'h 


'h 


'h 


,-0.08 
'h 



NOTE.- "The ID phase corresponds to our model with [3 = 
and o = 0. 6 The 2D phase (b) corresponds to our model with 
P = 1.1 and a = 0. 



simulations. Below we compare of the present work with the 
hydrodynamic simulation of S02 in Table 1. 

In the "ID" phase, the results of S02 can be well de- 
scribed by our model with /3 = and a — 0. Note that this 
"ID" phase corresponds to the evolutionary model with con- 
stant Ah (BC89). For «hs, the power law index is slightly 
(~ 10%,) different from our model (also BC89) and the re- 
sults of S02. In this case, P c oc l^ 1 and Phs = const are pre- 
dicted by this work and BC89, which coincides with the nu- 
merical results of S02 (see Fig. 6 (c) for P c and Phs in S02). 
In addition, our model can reproduce the constant pj (see 
Fig. 5(a) in S02). For comparisons, let us briefly comment 
on the self-similar model of expanding cocoons in which the 
growth of the cocoon head is included (e.g., Begelman 1996: 
hereafter B96). As already pointed out (e.g., Carvalho & 
O'Dea 2002), the self similar mofel of B96 cannot represent 
the behavor of the "ID" phase. The behavior of P c /Phs is 
also the intriguing issue. The decrease of P c /Phs with time 
is reported in Fig. 6 of S02. Using our model, this behavior is 
clearly explained by the decrease of the cocoon aspect ratio 
(see eq. (20)). 

The "2D" phase of S02 is well described by our model 
with (3 = 1.1 and a = 0. We adopt f3 = 1.1 to reproduce the 
P c evolution in Fig. 6 (c) of S02 because the other quantities 
shows much larger fluctuations in Fig. 6 of S02. The present 
model predicts the evolution of the hot spot pressure and 
mass density of the jet as Phs oc 1^ «hs c< Zj~ ' 56 and p s oc 
' 1 . These coincides with the average value of Phs, «hs, 
and pj (see Fig. 5 (a) and Fig 6 in S02). In the "2D" phase, 
the cross section of cocoon head grows as Ah oc l^' 1 unlike 
the "ID" phase (Ah = const). Thus, the velocity of hot spot 
decreases with Zh. Actually, the growth of the cross section 
area of the cocoon head can be seen in their simulations (Fig. 
3 in S02). In this phase, B96 also explains these results of 
S02. Moreover, the cocoon pressure is proportional to Phs 
in this phase of S02. From eq. (20), it can be understood 
with a constant 1Z. From above detailed comparison with 
"2D" relativistic hydrodynamic simulations, we found that 
the model represented in this paper can describe the flow 
and cocoon behavors seen in the "ID" and "2D" phases 
very well. It should be stressed that our analytic model is 
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Table 2. Comparisons with observations 





VHS 


Phs 


ms 


Pus 


Observations" 


const 


,-(1.3-1.7) 

'h 


jO.7-1.3 
'h 


,-(1.9-2.9) 
'h 


This work 6 


const 


,-(1.5-2.0) 


;0.75-l 


,-(1.5-2.0) 
'h 



NOTE.- a The results are adopted from Readhead et al. (1996a), 
Jeyakumar & Saikia (2000) and Pcrucho & Marti' (2002). 6 We set 
the slope index of the ambient density a = 1.5-2. 



more useful than numerical simulations when investigating 
a longer-term evolution of jets. The length of jets performed 
by numerical simulations of jets achieves at most the length 
order of 100 times of a jet beam size, while the spacial sizes 
of actual jets in AGNs are spread in six order of magnitude 
(i.e., from parsec to mega-parsec scale). 



3.2 Comparison with observations 

Based on a number of recent reports of indicating that the 
constant speed of hot spot advance (0.01 < «hs/c <0.1) 
(e.g., Readhead et al. 1996b; Carilli et al. 1991; Conway 
2002), we here examine the case of vhs = const. Obser- 
vationally, Phs and phs are inferred by using the mini- 
mum energy assumption and the neglecting the effect of 
thermal components (Readhead et al. 1996a; Jeyakumar & 
Saikia 2000; Perucho & Marti' 2002). From eq. (15), the 
relation of 2 — X(2 — 0.5a) = is required for the con- 
stant hot spot velocity. By eliminating the parameter X, 
our model reduce to the following forms; v c oc l h ^ 2 ^^ a 4 \ 
P c oc ^ /(a ~ 4) , Phs oc l~ a , p H s oc l~ a , r H s oc l^ /2 , pj oc l~ a , 
and TZ oc /-(«- 2 )/(«- 4 ). Here we used mean density pro- 
files obtained by a large number of sample clusters of galax- 
ies, which is pa(d) oc d~^ 1,5 2 - ) (e.g., Mulchaey & Zabludoff 
1998). 

We show the comparison with observational data for 
CSOs and FR II sources in Table 2. This indicates that 
our model well reproduce observational trends within the 
error bars. These agreements are likely to support "youth 
radio source scenario" basically. At the same time, the large 
dispersion of the observational data could tell us other pos- 
sibilities of evolutionary tracks of radio loud AGNs usually 
discussed. To explore it, it must be valuable to inquire into 
the origin of their large dispersion. Furthermore, we em- 
phasize that the deviation from the self similar evolution 
are frequently indicated by several authors (e.g., De Young 
1997, Gilbert et al. 2004). 

It is worth to show the reliability of the relation of the 
opening angle of hot spots derived by eq. (21), namely #hs = 
(vhs/c) 1 ^ 2 ^. For this aim, we adopt this equation to an 
archetypal radio galaxy Cygnus A. Using the values of vhs ~ 
0.01c (Carilli et al. 1991) and TZ = 0.6 (Wilson, Young, 
& Shopbell 2000) our model predicts #hs — 0.036, while 
the direct estimate of #j with rus = 2 kpc and z = 60 kpc 
indicate Ous — 0.033. Thus we can verify the reliability of 
the relation for the opening angle of hot spots and then we 
propose a new way of the estimation of vhs from obsevable 
two quantities TZ and #hs- It would be worth to compare 
with evaluations from the kinematic studies. 



4 CONCLUSIONS 

In the present work, we model a dynamical evolution of hot 
spots in radio loud AGNs. In this model, the unshocked flow 
satisfies the conservations of the mass, momentum, and ki- 
netic energy. We take account of the deceleration process of 
the jet by shocks, and the cocoon expansion which is iden- 
tified as the by-product of the exhausted flow. The model 
describes the evolution of physical quantities (vns, Phs, and 
Phs) in the hot spot in terms of lh- Below we summarize the 
main results based on this model. 

(i) We find that the ratio of P c /Phs is controlled by the 
aspect ratio of the cocoon TZ and slope index of the ambient 
medium a. If TZ remains to be constant during the jet prop- 
agation, the value P c /Phs is proportional to TZ 2 with the 
coefficient of order unity. This naturally explain the basic 
concept of the elongated cocoon with larger Phs than P c . 
Concerning P c , it is proportional to the bulk kinetic power 
of the jet in given p a . This is originated from our treatment 
of adiabatic injection of the dissipation energy of the jet into 
the cocoon. In addition, we suggest a new method to evalu- 
ate the velocity of hot spots from the aspect ratio of cocoon 
and the opening angle of hot spots. 

(ii) Our analytic model can well explain the results of 
2D co-evolution of jets and cocoons obtained by relativistic 
hydrodynamic simulations. This clearly guarantees the reli- 
ability of our model during the over-pressure cocoon phase. 
Since the dynamical length of jets obtained by numerical 
simulations is about a few 100 times of the jet beam size, our 
analytic model must be an useful tool to explore a longer- 
term dynamical evolution of jets than this scale. 

(iii) Our model prediction reasonably coincides with the 
recent observational trends of hot spots seen in CSO and 
FR II sources. Furthermore, we predict TZ oc l^ ' 2 °' and 
Ah oc li' 5 ~ 2 although little is done about systematic studies 
on these quantities. 

Lastly we should keep in mind that the present model do 
not take account of the details of (i) the absolute value of 
the mass density of the ambient medium, and (ii) radiative 
cooling effect which may be important for younger radio 
galaxies. In order to investigate whole story of evolutionary 
track of the radio loud AGNs, the study of above two points 
will be inevitably required. We plan investigate both of them 
in the near future. 
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